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I. INTRODUCTION 

Recent experimentsA*2iiM£ have shown that magnetic- 
field dependences of thermodynamic quantities in the 
vortex state of type-II superconductors provide unique 
information on the pairing symmetry and gap anisotropy. 
On the theoretical side, however, calculations of those 
quantities still remain a fairly difficult task to perform. 
The quasiclassical equations derived by Eilenberger— pro- 
vide a convenient starting point for this purpose. PeschX 
obtained a compact analytic solution to them based on 
the lowest Landau-level approximation for the pair po- 
tential. It has been used recently to discuss the field 
dependences of the thermal conductivity^ the density 
of states^ii2iii and thermodynamic quantities. 12 On the 
other hand, Kleini&±4 obtained a full numerical solution 
for the vortex-lattice state using a standard procedure 
to solve ordinary differential equations. This numeri- 
cal approach has been used extensively by Ichioka et 
a?i 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 to clarify the field dependences 
of the pair potential and the density of states for type-II 
superconductors with various energy-gap structures. It 
should be noted, however, that those numerical studies 
all adopted simplified model Fermi surfaces instead of 
complicated Fermi surfaces for real materials. Indeed, 
recent theoretical studiesifli 2 ^ 2 ^^ have clarified that de- 
tailed Fermi-surface structures are indispensable for the 
quantitative description of the vortex state. 

With these backgrounds, we here develop an alterna- 
tive method to solve the Eilenberger equations for the 
vortex-lattice states. A key point lies in expanding the 
pair potential and the quasiclassical / function in the ba- 
sis functions of the vortex-lattice states, thereby trans- 
forming the differential equations into algebraic equa- 
tions. This method has been powerful for (i) solving 
the Ginzburg-Landau equations^ and the Bogoliubov- 
de Gennes equations22*22iiffii2I and (ii) obtaining quan- 
titative agreements on the upper critical field H C 2 of 
Nb, NbSe2, and MgB2 with Fermi surfaces from first- 
principles electronic-structure calculations Thus, 



the method may be more advantageous for the cal- 
culations of thermodynamic quantities in finite mag- 
netic fields when realistic Fermi surfaces are used as in- 
puts. It will also be convenient for microscopically cal- 
culating responses of the vortex-lattice state to external 
perturbations^ such as vortex-lattice oscillations. 

This method is applied here to calculate magnetic-field 
dependences of the magnetization, the entropy, the Pauli 
paramagnetism, and the pair potential at various temper- 
atures for the two-dimensional s- and c? a .2_ y 2-wave super- 
conductors in the clean and dirty limits. These quantities 
have been obtained near H C 2 for the s-wave pairingiSiM 
Our purpose here is to clarify the overall field dependence 
of those quantities. 

This paper is organized as follows. Section II gives the 
formulation. Section III presents numerical results. Sec- 
tion IV summarizes the paper. We use k^ = \ throughout. 



II. FORMULATION 
A. Eilenberger equations 

We take the external magnetic field H along the z axis 
and express the vector potential ag&SLSSSiSLS 

A(r) = Bxy + A(r) . (1) 

Here B is the average flux density produced jointly by the 
external current and the internal supercurrent, and A is 
the spatially varying part of the magnetic field satisfying 
J VxA<ir = 0. We choose the gauge such that VA = 0. 

The Eilenberger equation for the even-parity pairing 
without Pauli paramagnetism is given by&S^ 

(e n +^{g) + \nvvty= (a0+A</)) 5 . (2a) 

Here e n = (2n+l)irT (n = Q, ±1, ±2, • • • ) is the Matsubara 
energy with T the temperature, r is the relaxation time 
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by nonmagnetic impurity scattering, and 
the Fermi-surface average: 



denotes Here a and a* are the boson operators: 



(g) = dS, 



ff(e n ,k F ,r) 
(27r)3JV(0)|v F | ' 

with dS*F an infinitesimal area on the Fermi surface, ./V(O) 
the density of states per spin and per unit volume at 
the Fermi energy in the normal state, and v F the Fermi 
velocity. The operator d in Eq. l(2"ajl is defined by 

d^V-i^- A, 

with <f>o = hc/2e the flux quantum^ A(r) is the pair po- 
tential, and 0(kp) specifies the gap anisotropy satisfying 
(</>(k F )) = 1. Finally, the quasiclassical Green's functions 
/ and g for e n >0 are connected by g—(l — fP) 1 ' 2 with 
/t(e n ,k F ,r) = /*(e n ,-k F ,r). 

Equation i|2a|) has to be solved simultaneously with the 
self-consistency equation for the pair potential and the 
Maxwell equation for A, which are given respectively by 



A(r) In 



71=0 



A(r) 



(</ , ( k F)/(£n,k F ,r)) 



(2b) 



-V 2 A(r) = -i 



167r 2 e7V(0)T 



X^( V F-9( £ "' k F> r )> ) ( 2c ) 



with T c q the transition temperature for r = oo. 

Finally, the free-energy functional corresponding to 
Eq. J5J is given bj&^i 



(X7 x A) 2 T 
<lr r\ } +jV(0)]A(r)l 2 ln^- 



-2ttTN(0) 



n=0 



|A(r) 



(J(e„,k F ,r)) 



(3) 



where F n is the free energy in the normal state and / is 
defined by 



1 = A7 + A/ t + 2e„( 5 -l) 



4r 



n iMzl fi / t VF-9/-/vF-av t (4) 



2r 



2(5+1) 



Indeed, functional differentiations of Eq. 10 with respect 
to /, A, and A lead to Eqs. (2ajl. (|2bjl . and QgJ, respec- 
tively. 
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with l c = 


= ^J$q/2ttB and c\c\ 



v x 

Vy - 2niBx/<&o 



, (6a) 



are defined by 



A : 



V2tt/ C/ - - . 

_ z — — (CiA x + lC2Ay) , 



V F+ = C 2 V Fx + lCiV Fy , 



(6b) 



(6c) 



respectively. The constants (ci,ca) can be fixed appro- 
priately to make the subsequent calculations efficient. A 
convenient choice 2 ^ is to impose the condition that the 
gradient term in the Ginzburg-Landau equation be ex- 
pressed in terms of a^a without using aa and aW, i.e., 
the pair potential near T c be described in terms of the 
lowest Landau level only. Alternatively, one may follow 
Graser et <il. " to change them at every temperature and 
magnetic field so as to make the free-energy within the 
lowest-Landau- level approximation smallest. 

Using Eq. (|6a|) . we can make up a set of basis functions 
to describe arbitrary vortex-lattice structures as2£ 



2ttL 



Cia 2 \fKV 



eX P 



n=-7V f /2+l 



iq y \ y 



x exp 



x exp 



na lx / , 2 

1— — \ V + l c q x 



n a 



l 2 c V ' 2 

c\c 2 fx — l 2 q y — na\-. 
~~2 



1 



-.H 



N 



x — l 2 q y — na\ x 



(7) 



Here N = 0, 1, 2, • • • denotes the Landau level, q is an 
arbitrary chosen magnetic Bloch vector characterizing 
the broken translational symmetry of the vortex lattice 
and specifying the core locations, and V is the volume 
of the system. The quantities a\ x , a\ v and a 2 are the 
components of the basic vectors ai and a2 in the xy 
plane, respectively, with a2 ||y and a\ x a2 — 2-nl 2 c , Af 2 de- 
notes the number of the flux quantum in the system, and 

Hn(x) = e x2 (—■£-) e ~ x2 is the Hermite polynomial. The 
basis functions are both orthonormal and complete, sat- 
isfying aipNq = VN'ipN-i<i and a^JVq = \/N + l^jv+iq- 



B. Operators and basis functions 



We first express the gradient operator in Eq. (|2a|l as 



•F + 



(a+A) - vp+i^+A*) 



(•5) 



C. Algebraic Eilenberger equations 

We now expand A, /, and A in the basis functions of 
the vortex lattice as 



A(r) = WJ2 A n VW r ) ; 



(8a) 



N=0 
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/(e n ,k P ,r) = VV^2 /Ar(e„,k F )'0Arq(r) , (8b) 



N=0 



l(r) = 



iKt 



(8c) 



where K is the reciprocal-lattice vector^ Substituting 
Eq. ijHJ) into Eq. (J2J) and using the orthogonality of the 
basis functions, Eq. (J2J) is transformed into a set of alge- 
braic equations for {/at}, {A^}, and {^4k} as 



and ff=(l — Z/^) 1 ^ 2 on those points by Eqs. JSJ) with re- 
stricting the summations to those satisfying N <N C and 
|K| <K C , where N c and K c are some cutoffs. Now, the 
integrations can be performed by the trapezoidal rule; its 
convergence is excellent for periodic functions. Also, the 
summation over n in Eq. (|9b|> is restricted in the actual 
calculations to those satisfying \e n \ < £ c - The conver- 
gence can be checked by increasing Ni nt , N c , K c , and 
e c . Finally, the vortex-lattice structure can be fixed by 
requiring that the free energy be minimum. 



E. Thermodynamic quantities 



J N + (3*VN+lf N+1 - PVNfN-i 

LJ ^ q Ltjy+h {f )g ~ <g)/ -p*A+pA*\ dr , (9a) 



Once A, /, and A are determined as above, we can 
calculate thermodynamic quantities of the vortex-lattice 
state. Specifically, the magnetization M due to super- 
current and the entropy S s are obtained by^ 



Aivln4^ = 2irT^ 



n=0 



— - (0(k F )/w(£™,k F )) 

£rj. 



IK 



167r 2 7V(0)T 
(Kl c B) 2 V 



with 



, (9b) 



/ (/3ff(£n,k F ,r))e- jK - r dr, 

71 = J 



(9c) 



-4ttM = / dr 

BV 



(VxA) 2 



.2^.VM.;7X( /V| -^, ( vH) ' r ). ,n,u 



n=0 



5 + 1 



S s — >S n 



N(0) 
T 



|rfr||A(r) 



-2ttT J2 [<*(£»> k F , r)) + 2e n (g-l)} I , (lib) 

n=0 ' 



Together with the equation to determine H C 2 derived 
recently^ the above coupled equations form a basis 
for efficient numerical calculations of the Eilenberger 
equations for vortex-lattice states with arbitrary Fermi- 
surface structures. 



respectively, where S n = 2tt 2 VN(0)T/3 is the entropy in 
the normal state and / is defined by Eq. Also, when it 
is much smaller than the diamagnetism by supercurrent, 
the magnetization M s p due to Pauli paramagnetism can 
be calculated by24 



M sP = M nP 



^ T ±fdr d M n 



V 



n=0 



de n 



VxA' 2 



B 



(11c) 

where M n p = 2VN(0)fi 2 i B with /zp, the Bohr magneton. 
The quantity 



dg_ 

de n 



2(1 -fP) 1 ' 2 Vdende. 



D. Numerical procedures 

For a given vortex-lattice structure specified by the 
basic vectors ai and a 2 in Eq. Q, the coupled equa- 
tion © may be solved iteratively in order of Eqs. (|9a[l . 
(|9bl) , and (|9"c)l by adopting a standard technique to solve 
nonlinear equations^ A convenient starting point is to 
put Af = S N0 A(T)^1-B/H c2 , f N = 0, and i K = 
on the right-hand of Eq. I|9a|) . where A(T) is the angle- 
averaged energy gap in zero field. In this connection, it 
may be worth noting that the tridiagonal matrix con- 
structed from the coefficients of /jy on the left-hand side 
of Eq. H9a(l can be inverted analytically. 2 ^ The Fermi- 
surface integrals (• • • ) can be performed as described in 
Sec. IV of Ref. 25. In contrast, integrations over r in 
Eqs. l|9a|) and (|9"c|l may be carried out as follows: At the 
beginning of each calculation, we prepare ipNqij) and 
e lKr at equally spaced N int x N int discrete points in a 
unit cell. We then construct A(r), /(e„,k F ,r), A(r), This equation can be solved similarly as Eq. (j2a| 



in Eq. (|llc|) may be obtained either by numerical differ- 
entiations or directly from the equation of differentiating 
Eq. ill' a I) with respect to e n : 



n d(g) 
2t de n 



f 



(13) 
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III. RESULTS 

We now present numerical results for two-dimensional 
systems with the cylindrical Fermi surface which is placed 
in the xy plane perpendicular to H. We have considered 
a couple of energy gaps: 



xlO 4 



0(k F ) 



1 



Fx ' 



s-wave 

dx 2 — y 2 ~ 



(14) 



Then it is convenient to set c\ = c 2 = 1 in Eq. © . There 
is another parameter in the system corresponding to the 
Ginzburg-Landau parameter n. We have fixed it by 



(fic/2e)A(0) 
^4^(0) h 2 v\ 



10 
7 



s-wave 
d 2 ,2_ J( 2-wave 



(15) 



with A(0) the angle-avera ged energy gap at T = 0. It 
follows from Eq. (46) of Ref . |33 that Eq. (|15fl corresponds 
to k — 49 and 40 for the s- and dj.2_ y 2-wave pairings 
in the clean limit, respectively, and k = 1300 for the s- 
wave pairing with r = 0.01fi./A(0). The vortex-lattice 
structure has been fixed as hexagonal (square) for the 
s-wave (c? a; 2_j / 2-wave) pairing so that finite contributions 
in the expansion (|8a|l come only from N = 0, 6, 12, • • ■ 
(N = 0, 4, 8, • • • ) Landau levels. 

Equation 10 for the above model has been solved with 
the procedure of Sec. IID over 0.05H c2 <B<H c2 . The 
value N cu t is chosen as 30 (16) for the s-wave {d x 2_ y 2- 
wave) pairing. On the other hand, we have set e c = 
20T C (50T C ) at T = 0.9T C (0.3T C ). These values have 
been enough to get the convergence. Using A, /, and 
A thus obtained, we have calculated the magnetization 
by supercurrent, the entropy, and the magnetization by 
Pauli Paramagnetism by Eqs. (|lla|) . (jllbfl . and i|llc[l . 
respectively. 

Figure ^ presents magnetic-field dependence of the 
magnetization by supercurrent for the s- and d x 2_ y 2-w&ve 
pairings in the clean limit at several temperatures. In 
both cases, the initial slope at B = H c2 gradually de- 
creases as the temperature is lowered, implying a mono- 
tonic increase of the Maki parameter— k 2 (T) as T^0. 
This feature of k 2 (T) has also been predicted in the case 
of the three-dimensional spherical Fermi surface with 
the s-wave pairing. 33 - 41 - 42 - 43 Unlike the three dimensional 
cnse^J^iiii however, the slope in these two-dimensional 
cases remains finite and does not approach even in the 
clean limit of T — > 0, in agreement with a previous cal- 
culation of The curves at low temperatures become 
more and more concave upward, thereby compensating 
the initial reduction of the magnetization. The tempera- 
ture variation is slightly larger for the d x 2_ y 2--wave pair- 
ing than the s-wave pairing. 

Figure [21 shows the expansion coefficients A at in Eq. 
(|8a|l as a function of B/H c2 for the s- and d x 2_ y 2- 
wave pairings in the clean limit at T = 0.3T C . Com- 
pared with the case near T c j2£ the mixing of higher Lan- 
dau levels develops from higher fields. However, the 
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FIG. 1: The magnetization by supercurrent as a function of 
B/H C 2 in the clean limit, (a) s-wave; (b) d x 2_ y 2-wa,ve. The 
temperatures are T/T c = 0.9, 0.7, 0.5, and 0.3 from top to 
bottom in the high-field region. 



contribution is still ~ 0.1 Ao even around B = 0.1H c2 . 
This fact implies that the Pesch approximation^ is excel- 
lent down to H ~ 0AH c2 for the two-dimensional cases 
with the isotropic Fermi surface. This may not be the 
case for systems with complicated Fermi surfaces, how- 
ever, as suggested by the fact that there is already an 
amount of higher-Landau-levcl contributions at H c2 in 
those cases The Pesch approximation may be im- 

proved to some degree by the procedure of Graser et alm& 
to change c\ and c 2 in Eq. (|6all at every temperature and 
magnetic field so that the free-energy is smallest. 

Figures and ^ plot the entropy S s and the magne- 
tization M s p due to Pauli paramagnetism, respectively, 
as a function of B / H c2 for the s- and d x 2_ y 2--w&ve pair- 
ings in the clean limit at several temperatures. To see 



TABLE I: The exponent as of Eq. I|16a^ for the s- and d x 2_ y 2- 
wave pairings in the clean limit calculated by the best fit to 
the numerical data of Fig. 3. 



T/Tc 


0.9 0.7 0.5 0.3 


as (s-wave) 


0.73 0.69 0.66 0.63 


as (d-wave) 


0.72 0.66 0.59 0.52 
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FIG. 2: The expansion coefficients A at in Eq. 18aH as a func- 
tion of B/H C 2 in the clean limit at T = 0.3T C : (a) s-wave 
pairing where the curves are Ao, A6, and A12 from top to 
bottom; (b) d x 2_ y 2-wa,ve pairing where the curves are Ao, 
A4, and Ag from top to bottom near H C 2- 




m 2 



(b) 

FIG. 3: The entropy S s as a function of B/H C 2- (a) s-wave; 
(b) dj,2_ H 2-wave. The temperatures are T/T c = 0.9, 0.7, 0.5, 
and 0.3 from bottom to top. They are normalized to vary 
from 1 at B = H c2 to at B = 0. 



the field dependence clearly, the entropy is normalized 
by using S s o = S S (B = 0) and S n = S S (B = H c2 ) as 
(S s -S s0 )/(S n -S s0 ); it varies from 1 to for H c2 >B>0. 
The same normalization is adopted for M s p. All curves 
deviate upwards from the linear behavior oc B / H c2 and 
become more and more convex upward as T ^ 0. This 
tendency is more conspicuous for the d 3 .2_ y 2-wave pairing 
due to the residual low-energy density of states. To see 
the behavior quantitatively, we have fitted our numerical 
data by the formulas: 



B Y s 



(16a) 



TABLE II: The exponent a x of Eq. l(l6b> for the s- and 
d x 2_ y 2 -wave pairings in the clean limit calculated by the best 
fit to the numerical data of Fig. 4. 



T/T c 


0.9 0.7 0.5 0.3 


a x (s-wave) 


0.71 0.67 0.63 0.63 


a x (d-wave) 


0.70 0.63 0.56 0.49 



M s p-M sP o 
M n — M s 



sPO 



B 

H C 2 



(16b) 



To confirm the numerical results, we first estimated as 
and a x for the s-wave pairing by using only the data of 
0.85i? C 2 <B< Q.95H c2 . Although not presented here, the 
procedure excellently reproduced the values of a previous 
calculation near 7J c2 f24 as they should. 

Table I shows the exponent as obtained by the best 
fit to the data of 0.05ff c2 <B< 0.95H c2 in Fig. 3. The 
value of the d x 2_ y 2 pairing at T/T c — 0.9 is almost the 
same as the corresponding s-wave result. As the temper- 
ature is decreased, however, the d x 2_ y 2-w&ve exponent 
decreases more rapidly so that the curve in Fig. 3(b) be- 
comes more convex upward; this is due to the residual 
low-energy density of states of the d x 2_ y 2-w&ve pairing. 
Table II presents another exponent a x obtained by the 
best fit to the data of 0.05H c2 <B< 0.95H c2 in Fig. 4. 
Each value is fairly close to the corresponding one for 
as, as may be expected from the fact that both quan- 
tities probe the zero-energy density of states. In this 
context, Ichioka et alH calculated the field dependence 
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FIG. 5: Magnetization by supercurrent as a function of B / H C 2 
for the s-wave pairing with r = ft/A(0). The temperatures are 
T/T c = 0.9, 0.7, 0.5, and 0.3 from top to bottom. 



1 - 8 " ^^<^ 
%" 

5 - 0.6- y^^y 

$ /// T/T=0.3 

3 0.4- — T/T=0.5 - 

\ /V/ -T/T=0.7 

g 0.2- byy t/t=o.9 ■ 

°'8.0 0.2 0.4 0.6 0.8 1.0 

B/H 2 

(b) 

FIG. 4: The magnetization M a p by Pauli paramagnetism as 
a function of B/H C 2- (a) s-wave; (b) d x 2_ y 2-wa,ve. The tem- 
peratures are T/T c — 0.9, 0.7, 0.5 and 0.3 from bottom to 
top. They are normalized to vary from 1 at B — H C 2 to at 
B = 0. 



of the zero-energy density of states iV(0) at T/T c = 0.5 
to find N(0) cx (B/H c2 ) - 67 and N(0) cx {B/H c2 f A1 for 
the s- and d x 2_ y 2-wave pairings, respectively. Our esti- 
mates for ots(T— +0) and a x (T— >0) are somewhat smaller 
(larger) for the s-wave (d a .2_ l ,2-wave) pairing. 

Experiments on the T-linear specific-heat coeffi- 
cient 7s(-B) have been performed for clean V3Si, 44 
NbSe 2) 45 ' 46 i 47 ' 48 and CeRu 2 i2 The quantity j s (B)/~/ n 
coincides for T -> with (S s - S s0 )/(S n - S s0 ) of Fig. 
3. Those data all show marked upward deviations from 
the linear behavior "f n B/H C 2, indicating that it is a com- 
mon feature among clean superconductors irrespective of 
the energy-gap symmetry. Sonier et a/4£ thereby ex- 
tracted the exponent 0.66 for the field dependence of 
7s(-B) as T— >0, in good agreement with the result 0.67 
by Ichioka et alH for the clean two-dimensional s-wave 
model. However, a more recent experiment by Hanaguri 
et alm£- reported a different exponent 0.5. It should also 
be noted that NbSe 2 has three kinds of Fermi surfaces 
and one of them is quite different in structure from the 
cylinder. There also exists a recent experiment which 
indicates existence of different superconducting energy 



gaps on different Fermi surfaces £ Hence the agreement 
between the experiment by Sonier et tdM* and the the- 
ory by Ichioka et alw^ might be an artifact and should be 
confirmed by more detailed experiments as well as theo- 
ries incorporating both Fermi-surface and gap structures. 
In this context, it is worth noting that no detailed exper- 
iments have been performed on the field dependence of 
7s(-B) even for the classic type-II superconductors V and 
Nb, although early experiments£2*£i suggest similar up- 
ward deviations from the linear behavior j n B/H C 2- 

We next focus on the s-wave pairing in the dirty limit. 
Figure [S] shows the magnetization for r = 0.01/i/A(0) as 
a function of B / H C 2- Compared with the clean- limit re- 
sults of Fig. 1(a), we observe an extended linearity down 
to B/H c2 ~ 0.2 irrespective of the temperature. The de- 
crease of the initial slope for T — > is as expected from 
the temperature dependence of the Maki parameter 
By scaling this change of the initial slope, all the curves 
almost fall onto a single curve. This is a marked feature 
in the dirty limit which is absent in the clean-limit result 
of Fig. 1(a). 

Figure [fj] shows the field dependences of S s and M s p 
for t = 0.01fi./A(0) at various temperatures. Table III 
presents the corresponding exponents a s and a x obtained 
from the data of 0.5iJ c2 < B < 0.95i? c2 ; unlike the clean- 
limit case, it has been impossible to fit the whole region 
by a single exponent, especially at intermediate temper- 
atures, as may be realized from Fig. Compared with 
Figs. 3(a) and 4(a), the curves are more monotonic with 
the almost linear behavior oc B / H C 2- Looking at the tem- 
perature dependence more closely, however, we observe 
a change from a convex-upward behavior at high tem- 



TABLE III: The exponents as and a x for the s-wave pairing 
with t = 0.01fi/A(0) calculated by the best fit to the numerical 
data of 0.5J?c2 <B< 0.95H c2 . 



T/T c 


0.9 0.7 0.5 0.3 


as 


0.84 0.88 0.92 1.10 




0.84 0.89 0.99 1.15 
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FIG. 6: (a) The entropy Ss and (b) the magnetization M s p 
by Pauli paramagnetism as a function of B/H C 2 for the s- 
wave pairing with r = 0.01fi/A(0). The temperatures are 
T/T c = 0.9, 0.7, 0.5 and 0.3 from top to bottom. 



peratures to a convex- downward behavior at low tem- 
peratures, in agreement with a previous calculation near 
H c2 ~ This feature also appears in the field dependence 
of the zero-energy density of states as calculated recently 
by Miranovic et aim! The convex-downward behavior at 
T = 0.3T C may become more pronounced at lower tem- 
peratures to be observable experimentally. 



main analytic formulas are given in Sec. IIC together with 
the numerical procedure to solve them in Sec. IID. This 
method, which directly extends the H C 2 equation2Si2& 
to lower fields, has a potential applicability to systems 
with complicated Fermi surfaces and/or gap structures to 
carry out detailed calculations on the field dependences 
of thermodynamic quantities for various type-II super- 
conductors. 

Using it, we have calculated the field dependences of 
the magnetization by supercurrent, the mixing of higher 
Landau levels in the pair potential, the entropy, and 
the Pauli paramagnetism for the two-dimensional s- and 
d x 2_ y 2-wa,ve pairings in the clean and dirty limits at var- 
ious temperatures. Previous results near H C 2 for the s- 
wave pairing22i2i have been reproduced adequately and 
extended to lower fields to clarify the overall field depen- 
dences. The differences between the s- and d x 2_ y 2-w&ve 
pairings are quite small at high temperatures but de- 
velop gradually as the temperature is lowered, reflecting 
the marked difference in the low-energy density of states 
between the two cases. The field dependences of the en- 
tropy and Pauli paramagnetism in the clean limit at low 
temperatures present convex-upward behaviors for both 
pairings. In contrast, the curves of the s-wave pairing in 
the dirty limit are more monotonic and fairly close to the 
linear behavior, but also acquire downward curvature at 
low temperatures. As for the magnetization by super- 
current, there is a wide region of linear field dependence 
from H C 2 both at high temperatures and in the dirty 
limit. The region shrinks in the clean limit as the tem- 
perature is lowered, and the curve acquires pronounced 
upward curvature. It is also found that the mixing of 
higher Landau levels in the pair potential is small for 
B > 0.1H c2 but develops rapidly as the field is further 
decreased. 
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